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We demonstrate that period-doubled discrete breathers appear from the anti-continuum limit of 
the driven Peyrard-Bishop-Dauxois model of DNA. These novel breathers result from a stability 
overlap between sub-harmonic solutions of the driven Morse oscillator. Sub-harmonic breathers 
exist whenever a stability overlap is present within the Feigenbaum cascade to chaos and therefore 
an entire cascade of such breathers exists. This phenomenon is present in any driven lattice where 
the on-site potential admits sub-harmonic solutions. In DNA these breathers may have ramifications 
for cellular gene expression. 
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Introduction. — Discrete breathers are spatially local- 
ized, temporally periodic excitations in nonlinear lattices 
[H . While discrete breathers share many traits with soli- 
tons they stand out because of their localization brought 
about by a delicate sensitivity to lattice discreteness. Dis- 
crete breathers, have been ubiquitously studied in a wide 
variety of physical systems and have been the subject 
of intense theoretical and numerical scrutiny [H-lHl- Al- 
though, dissipation and driving are typically key exper- 
imental features, most of the theoretical and numerical 
studies have excluded such effects. Only quite recently 
have the effects of dissipation and driving been explicitly 
considered in numerical studies of discrete breathers 0. 

Biomolecules represent a striking example of systems 
where the localization of discrete breathers has long 
been considered important, and where dissipation and 
driving cannot be ignored. DNA represents a spe- 
cific biomolecular system where localization of energy 
in terms of strand-separation dynamics is emerging as 
a governing regulatory factor [7|. The Peyrard-Bishop- 
Dauxois (PBD) model of double-stranded DNA d |9| is 
arguably the most successful model for describing this lo- 
cal pairing/unpairing (breathing) dynamics, since it re- 
produces a wide varietyof experiments related to strand- 
separation dynamics [10|. Of particular importance is 
the role this model has played in demonstrating strong 
correlations between regulatory activity, such as protein 
binding and transcription, and the equilibrium propen- 
sity of double-stranded DNA for local strand separation 
ll-l3|- Recently, this model was augmented to include 
a monochromatic drive in the terahertz frequency range, 
and was suggested to represent a simplified model for 
DNA dynamics in the presence of terahertz radiation 
[l5| . In this recent work [l^ the existence of discrete 
breathers oscillating at half the frequency of the drive 
was numerically observed, and it was argued that the 
existence and spontaneous generation of such breathing 
states may have significant ramifications for cellular gene 
expression: This argument has received some experimen- 
tal support Til 17|. 



Here, we investigate these period-doubled breather ex- 
citations further, and show that they appear naturally 
from the anti-continuum limit [l8| in systems of nonlin- 



ear oscillators where multi-stability occurs between var- 
ious stages of the Feigenbaum period doubling cascade. 
Although, our study focuses on the PDB model the ob- 
served phenomena are rather general and not limited to 
this specific system of nonlinear oscillators. 

DNA model. — We consider the PBD model in the fol- 
lowing form [l5l |: 

myn = - V'iyn) - W(yn+\,yn) " VF' (?/„ , y„_ i ) 

- m7y„ -f Focosjli, (1) 

where the Morse potential 

lJ{yn) = D [exp {-ayn) - if 

represents the hydrogen bonding of the complementary 
bases. Similarly, 

W{yn,yn-i) - ^ (l + pe-^<^y-+y"-^^) (2/„ - y„_i)' 

represents the stacking energy between consecutive base 
pairs. The parameters D, a, and k in principle de- 
pend on the type of the base pair (A-T or G-C), but 
for simplicity we will consider homogeneous DNA in this 
study. The term mjyn is the drag caused by the solvent, 
while Fq cos fit is the (terahertz) drive [l^. In this sys- 
tem, Eq. ([1}, the linear resonance frequency is given as 
nl = 2Da^/m. 

In order to understand how period-doubled breathers 
oscillating at half the frequency, of the drive can 

appear at the anti-continuum limit, we first study the 
nonlinear response of a single {k = 0) oscillator. For 
periodic driving the single oscillator supports periodic 
solutions with the period (T = 27r/f2) of the drive, and, 
due to the softness of the Morse potential, this oscillator 
also supports solutions with period nT, where n is any 
positive integer. These solutions are most simply studied 
at zero dissipation (7 = 0); we will then examine the 
modifications resulting when 7 7^ 0. A Newton method 
was used to track the solutions with periods T and 
2T. The Newton procedure was initiated for a very small 
driving amplitude {Fq ~ 0). The response of the system, 
namely the amplitude of the period T and 2T solutions, 
was followed as li^ol was increased. 
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FIG. 1; (Color online) Nonlinear response manifold of a sin- 
gle oscillator driven at Q =7.0 rad/ps (< Qo)- The NLRM 
for the period-T solutions (solid black line) and the period- 
doubled solutions (dashed red line) are shown together with 
parts of the NLRM for the period-4r (solid blue line seg- 
ments connected to BP4) solutions. The stability of the vari- 
ous solutions is indicated by the letters placed in the vicinity 
of the lines ('s' for stable and 'u' for unstable). Important 
turning points {TP), bifurcation points (BP), and crossing 
points (CP) are shown by circles. Finally, the grey areas rep- 
resent the chaotic regions that result from the Feigenbaum 
period-doubling cascades. The inset shows the behavior near 
CPi in the presence of dissipation. 7 = (dashed red line), 
7 = 0.001ps-\ 7 = 0.01ps-\ and 7 = Clps'^ 
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Soft period- doubled breathers: Q < CIq. — In the above 
fashion we constructed the nonlinear response manifolds 
(NLRM) ^] depicted for r2 = 7 rad/ps in Fig. [H The 
NLRM for the oscillations with period T is shown by 
the solid black line, while the dashed red line represents 
the NLRM for the period-doubled {2T) solutions. Fi- 
nally, the blue lines segments (connected to BP4) rep- 
resent parts of the NLRM for the period 4T solutions. 
The stability of the solutions can be obtained by Floquet 
stability analysis and is indicated by the letters 's' 
(stable) and 'u' (unstable) placed adjacent to the respec- 
tive NLRMs. The stability of a given solution changes 
at the turning points (TP), bifurcation points (BP), or 
at the crossing points (CP) of the NLRM. Since, the 
Morse potential is asymmetric with respect to y = 0, the 
NLRMs lack symmetry with respect to a 7r-phase shift 
(±i^o) of the drive. However, an equivalence between the 
NLRM branches for Fq and —Fq remains, as is indicated 
by our labeling of the TPs, BPs, and CPs. 

From Fig. [1] it can be noticed that, for a region of 
small positive force amplitudes, Fq, the period-T solu- 
tion is stable up to TPi (Fo ~ ±14 pN) and so is the 
period-doubled solution between CP2 {Fq — 0) and BP3 
(Fo ~ 21.5 pN). This means that in the region between 
< Fq < 14 pN both solutions are stable. A similar 
stability overlap can be observed for the regions Fq = 
to TPi (period-T, Fq < 0) and BP4 to BPi (period- 



FIG. 2: Illustration of period-doubled breathers. Shown are 
the central site yo of the breather (solid red lines) together 
with the closest neighbor sites j/i (dashed blue lines), and 
y2 (dashed-dotted green lines). For clarity, panel a) shows 
lOyi and lOj/2, respectively. Panel a) shows the period- 
doubled breather corresponding to the stability overlap shown 
in Fig. [1] (n = 7 rad/ps and 7 = 0) for Fo = 7.06 pN and 
k = 0.02eVA~'^. Similarly, panel b) shows a period-doubled 
breather for = 5.35 rad/ps, 7 = lps~^, Fo = 115.6 pN, and 
k = 0.0065eVA~'^. The insets shows the full spatio-temporal 
evolution of the period-doubled breathers. 



doubled) when fl < 5.5 rad/ps (not shown in Fig. [TJ. 
The ramifications of these stability overlaps for the con- 
struction of discrete breathers is now clear. In the anti- 
continuum limit (fc = 0) of driven oscillators, we can 
construct solutions where — 1 of the oscillators perform 
stable period-T oscillations, while the remaining oscilla- 
tor resides stably in the period-doubled state. Perform- 
ing standard numerical continuation [2l[ of this state to 
finite values of the coupling, fc, will retain the stabil- 
ity of the solution and lead to a stable period-doubled 
discrete breather solution. Examples of such period- 
doubled breathers arising in each of the two described 
stability overlap regions are shown in Fig. [5] Panel a) of 
Fig. [5] shows the period-doubled breather in the stabil- 
ity overlap region for Fq > 0. It is clearly seen that the 
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FIG. 3: Nonlinear response manifold of a single oscillator 
driven at =13.194 rad/ps (> f2o). The NLRM for the 
period-T solution (solid black line) and the period-doubled 
solutions (dashedred line) are shown together with parts of 
the NLRM for the period-4T (blue line segments) solutions. 
The stability of the solutions is indicated by the letters placed 
in the vicinity of the lines ('s' for stable and 'u' for unstable). 
Important turning points, bifurcation points, and crossing 
points are shown by circles. Finally, the grey areas represent 
the chaotic regions that result from the Feigenbaum period- 
doubling cascades. A stability overlap can be observed in the 
vicinity of Fq — 0. The inset shows a magnified version of the 
area between the chaotic regions. 



breather, which extends over 5 lattice sites, performs os- 
cillations at exactly half the frequency of the background. 
The background is frequency-locked to the drive at fre- 
quency O. Similarly, panel b) of Fig. [5] shows the period 
doubled breather that exists in the stability overlap re- 
gion for Fq < 0. The breather is more localized and 
its dynamical behavior is more complex. This breather 
roughly follows the background period-T oscillations and 
the period-doubling occurs as a modulation of the breath- 
ing amplitude. For weak driving the dynamics must be 
very similar to the undriven case, as is seem in panel a). 
At stronger driving the dynamics is more dominated by 
the drive and become more intricate, as in the case of 
panel b). 

It should be noted that, in the absence of a stability 
overlap, unstable period-doubled breathers can of course 
still be constructed. An interesting example of this oc- 
curs for Fq < 23.3 nN where Fig. [T]shows that the period- 
doubled solution is stable while the period-T solution is 
unstable. Close to BP2 the period-T solution is very 
weakly unstable (the Floquet exponent is ~ 1.04), so that 
a weakly unstable 'dark' breather can be constructed in 
this region by positioning N-1 of the oscillators in the 
period-doubled states while the remaining oscillator is in 
the period-T state. We have found that numerical con- 
tinuation of this state produces a weakly unstable 'dark' 
breather that can be sustained for tens of periods. 



The scenario described above generally holds also when 
dissipation is taken into account. The modifications to 
the NLRM caused by dissipation is thoroughly discussed 
for a different system of nonlinear equations in Ref. . 
The main effect of dissipation occurs at the crossing 
points CPi and CP2- For 7 = the period-T solutions 
are in phase with the drive when the response amplitude 
Ho is below CPi, but it is in anti-phase (7r-phase shifted) 
when j/o is above CPi. Similarly, the period-doubled so- 
lutions experience a 7r-phase shift at CP2. The introduc- 
tion of dissipation creates an additional phase shift be- 
tween the drive and the response everywhere along the 
NLRM. This additional phase shift is largest close to the 
crossing points and cause the solution to disappear close 
to the crossing points (small \Fo\). This effect is illus- 
trated in the inset of Fig. [T]for the CP2 crossing point. 
It is clear that for weak dissipation the described stability 
scenario persists. However, as the dissipation increases 
the stability overlap region is eventually lost and with it 
the stable period-doubled breathers. 

Hard period- doubled breathers: > flo. — Figure 
[3] shows the NLRM when the frequency of the drive 
(ri = 13.194 rad/ps) is larger than the linear resonance 
frequency of the on-site potential. In this case only a 
single, but larger, stability overlap between the period- 
T solution and the period-doubled solution exists in the 
vicinity of Fq = 0. Again, the period-doubled breathers 
can be generated in this overlap region. It is noteworthy 
that for soft potentials, such as the Morse potential, or- 
dinary period-T breathers cannot exists above the linear 
f2o resonance frequency: This frequency region can only 
be accessed for hard potentials. 

Two examples of the period-doubled breathers for = 
13.194 rad/ps are given in Fig. |4l Panel a) shows a period 
doubled-breather in a dissipation-less (7 = 0) system, 
and it is clear that the breather sites oscillate at half the 
frequency of the background, which is frequency locked 
to the drive. Similarly, panel b) shows a period-doubled 
breather at strong dissipation 7 = Ips^^. The phase 
shift introduced by the dissipation is apparent. 

Discussion and conclusions. — We have demonstrated 
the existence of a new kind of sub-harmonic discrete 
breathers in driven lattices of coupled nonlinear oscil- 
lators. The phenomenon was illustrated within the 
driven PBD model, which consists of a lattice of coupled 
Morse oscillators. Specifically, we have shown that these 
breathers appear naturally from the anti-continuum limit 
in a fashion similar to familiar discrete breathers. We 
show that the sub-harmonic breathers arise from a stabil- 
ity overlap between sub-harmonic solutions of the single 
oscillator. This new kind of breathers exist even when 
the frequency of the drive is above the linear resonance 
frequency of the soft on-site potential. Ordinary breather 
solutions do not exist in such cases for soft potentials. 

Although, we focused here on the stability overlap 
between the frequency locked solution and the period- 
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double-stranded DNA, and thereby could lead to modifi- 
cation in cellular gene expression by affecting transcrip- 
tion and other processes. Such effects have recently been 
experimentally observed [l7j . Future extensions of 
our approach will include effects of gene sequence and 
other disorder, as well as temperature. 

This research was carried out under the auspices of the 
National Nuclear Security Administration of the U.S. De- 
partment of Energy at Los Alamos National Laboratory 
under Contract No. DE-AC52-06NA25396. 
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FIG. 4: Illustration of period-doubled breathers. Shown are 
the central site yo of the breather (solid red lines) together 
with the closest neighbor sites yi (dashed blue lines), and 
2/2 (dashed-dotted green lines). Panel a) shows the period- 
doubled breather corresponding to the stability overlap shown 
in Fig. [U (r2 = 13.194 rad/ps and 7 = 0) for Fo = 160 pN 
and k = 0.025eVA~^. Panel b) demonstrates the effects of 
dissipation = 13.194 rad/ps, 7 = lps~^, Fq = 209 pN, and 
k — 0.025eVA~^. The insets show the full spatio-temporal 
evolution of the period-doubled breathers. 



doubled solution, it is clear that in the presence of a sta- 
bility overlap the phenomena can occur between any of 
the sub-harmonic solutions in the Feigenbaum cascade. 
We have explicitly verified this by constructing stable 
periodic-nT (n = 3,4,5, and 6 ) breathers for small fo 
and fc, in the case where > Qq. In this fashion a 
cascade of sub-harmonic discrete breathers exists. The 
phenomenon is not specific to the system we have stud- 
ied: it is present in any driven lattice of nonlinear oscil- 
lators where the on-site potential admits sub-harmonic 
solutions. The sub-harmonic breathers also exist in the 
presence of dissipation, it is however important to note 
that in this case the cascade is truncated. 

The existence of the sub-harmonic breathers in DNA 
is potentially very important because this suggests that 
electromagnetic terahertz radiation can lead to the spon- 
taneous generation of stable local strand separation in 
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